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We compute the axial and pseudoscalar form factors of the nucleon in the Dyson-Schwinger ap- 
proach. To this end, we solve a covariant three-body Faddeev equation for the nucleon wave function 
and determine the matrix elements of the axialvector and pseudoscalar isotriplet currents. Our only 
input is a well-established and phenomenologically successful ansatz for the nonperturbative quark- 
gluon interaction. As a consequence of the axial Ward-Takahashi identity that is respected at the 
quark level, the Goldberger-Treiman relation is reproduced for all current-quark masses. We discuss 
the timelike pole structure of the quark- ant iquark vertices that enters the nucleon matrix elements 
and determines the momentum dependence of the form factors. Our result for the axial charge 
underestimates the experimental value by 20 — 25% which might be a signal of missing pion-cloud 
contributions. The axial and pseudoscalar form factors agree with phenomenological and lattice 
data in the momentum range above ~ 1 . . . 2 GeV^. 
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I. INTRODUCTION 

The nucleon's axial and pseudoscalar form factors 
are of fundamental significance for the properties of 
the nucleon that are probed in weak interaction pro- 
cesses. Their momentum dependence can be experimen- 
tally tested by (anti)neutino scattering off nucleons or 
nuclei, charged pion electroproduction and muon capture 
processes; see [1-3] for reviews. Both form factors are ex- 
perimentally hard to extract and therefore considerably 
less well known than their electromagnetic counterparts. 
Precisely measured is only the low-momentum limit 
of the axial form factor which is determined from neu- 
tron /?-decay. Planned experiments at major facilities are 
expected to change this situation in the near future. 

The theoretical calculation of the nucleon's axial 
and pseudoscalar form factors requires genuinely non- 
perturbative methods. Chiral perturbation theory has 
been successful in this respect [1, 4, 5] although it is gen- 
erally limited to the region of low momentum transfer. 
Recent studies in lattice gauge theory are getting closer 
to the physical pion mass region [6-8] but finite-volume 
effects become increasingly important. Another non- 
perturbative approach is the one via functional meth- 
ods, in particular Dyson-Schwinger equations (DSEs) 
and Faddeev equations. The basic idea here is to start 
at the level of the Green functions for quarks and gluons 
and construct hadronic bound states from correspond- 
ing Bethe-Salpeter equations (BSEs) and Faddeev equa- 
tions [9-11]. Those bound states are subsequently probed 
by the respective currents to yield the form factors of 
interest. One of the important advantages of this ap- 
proach is its direct access to the substructure of mesons 
and baryons at all momentum scales and values of the 
current-quark mass. In principle such a setup is ideal to 
identify mechanisms that are responsible for the wealth 
of experimental phenomena associated with the study of 
hadrons. 



The study of axial and pseudoscalar form factors in 
the functional approach has so far been limited to an 
approximation where the nucleon is treated as a bound 
object of a quark and a diquark that interact via quark 
exchange [12, 13]. The entire gluonic substructure ap- 
pears here only implicitly within the dressing of quark 
and diquark propagators as well as diquark vertex func- 
tions. There are several conceptual issues that compli- 
cate the treatment of form factors in the quark-diquark 
model. First, the requirement of current conservation in- 
duces the appearance of intricate 'seagull' diagrams [14]. 
Such terms have been taken into account for electromag- 
netic form factors, but their implementation in the case 
of axial form factors has not yet been possible for tech- 
nical reasons [13]. Second, to comply with chiral Ward 
identities, a current-conserving quark-diquark model re- 
quires vector diquarks in addition to the usual scalar and 
axialvector diquark degrees of freedom [15] . Such an elab- 
orate treatment of the quark-diquark model has not yet 
been performed. 

The situation is somewhat different when the nucleon 
is treated as a genuine three-body problem. The re- 
sulting Faddeev equation in rainbow-ladder truncation 
has been solved only recently for the nucleon and A 
masses [16, 17], and the corresponding nucleon electro- 
magnetic form factors have been calculated in Ref. [18]. 
Here the currents couple to the fully dressed quark prop- 
agators only and thereby provide a clean access to the 
underlying quark and gluon substructure, without the 
additional complications of the quark-diquark model. In 
the present work we dwell on this advantage and deter- 
mine the nucleon's axial and pseudoscalar form factors 
in the three-body approach. 

The paper is organized as follows. In Sec. II we sum- 
marize the details of the three-body Faddeev framework 
and specify the effective quark-gluon coupling that serves 
as a model input in our calculation. We furthermore dis- 
cuss the structure of the axial and pseudoscalar vertices 
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that are central for our calculation of the respective form 
factors. Here we made additional progress in determining 
these vertices explicitly from their inhomogenuous Bethe- 
Salpeter equations instead of the model ansatze that were 
used in previous works [13]. In Sec. Ill we briefly dis- 
cuss the various interrelatons between the axial, pseu- 
doscalar and pion-nucleon form factors and the resulting 
Goldberger-Treiman relation. Our results are presented 
in Sec. IV and we conclude in Sec. V. Throughout this 
paper we work in Euclidean momentum space and flavor- 
SU{2) in the isospin-symmetric limit m„ — m^. 



II. FADDEEV FRAMEWORK 

Our goal is to compute the matrix elements of the ax- 
ialvector and pseudoscalar isotriplet currents 



q{x) 757'' y q{x) 



and q{x) 175 y q{x) (1) 



in the nucleon. The ak represent the three Pauli matrices 
in SU{2). The nucleon matrix elements obtained from 
Eq. (1) are expressed in terms of three Lorentz- invariant 
form factors that depend on the momentum transfer Q^. 
The axial form factor Ga{Q^) and induced pseudoscalar 
form factor Gp{Q^) parametrize the axialvector current, 
whereas the pseudoscalar form factor G^{Q^) determines 
the pseudoscalar current. We will discuss the respective 
decomposition in detail in Section III. 

The computation of these form factors in the Dyson- 
Schwinger framework necessitates a quark-level picture 
of the nucleon current matrix elements. In the follow- 
ing subsections we will outline a description where all 
ingredients of the current diagrams are computed self- 
consistently once the quark-(anti-)quark interaction is 
specified. The nucleon's bound-state amplitude as so- 
lution of the covariant Faddeev equation is discussed in 
Section II A and the ansatz for the kernel and corre- 
sponding equation for the dressed quark propagator in 
Section II B. The construction of the nucleon's current 
matrix elements is described in Section II C, and in Sec- 
tion II D wc detail the properties of the axialvector and 
pseudoscalar qq vertices that enter the form factor dia- 
grams. 



A. Covariant Faddeev equation 

The starting point of the covariant bound-state ap- 
proach is the three-quark Green function G, as well as 
its three-quark connected and amputated counterpart, 
the scattering matrix T that is defined via 



G — Go + Go T Go . 



(2) 



Here, Go = S^S^S denotes the disconnected contribu- 
tion to G, i.e., the product of three dressed quark prop- 
agators. In order to keep the discussion transparent, we 



will frequently resort to a symbolic notation where Dirac- 
Lorentz, color and flavor indices as well as momentum de- 
pendencies are suppressed, and the products in Eq. (2) 
and related equations are understood as four-momentum 
integrations over all internal loop momenta. The Green 
function G satisfies a scattering equation, which is the 
nonperturbatively resummed Dyson series 



G = Go + Goi^G ^ G- 



Gn 



(3) 



and the equivalent relation for the T— matrix reads 

T = K + KG0T <^ T-^=K-^-Gq. (4) 

The input in both equations is the three-quark kernel K 
that will be detailed below. 

Baryons in QCD appear as poles in the scattering ma- 
trix T. Such a pole defines a baryon bound state on its 
mass shell = — M^, where P is the total baryon mo- 
mentum and M is its mass. The scattering matrix at the 
pole assumes the form 



P2 + Af 2 



(5) 



and thereby defines the baryon's covariant bound-state 
amplitude 4*, with 5* being its charge conjugate. At 
the pole, Eq. (4) reduces to a homogeneous bound-state 
equation for the amplitude 



5- = A' Go * , 



K = K 



m 



(a) 



(6) 



The three-body kernel K contains a three-quark irre- 
ducible contribution K^^] and the sum of permuted two- 
quark kernels -ft^(a)j where the subscript a denotes the 
respective spectator quark. 

The first step in constructing a feasible ansatz for K is 
to disregard irreducible three-quark interactions such as, 
in the simplest case, those originating from a three-gluon 
vertex that connects all three quark lines. Omitting K^^^] 
in Eq. (6) yields the covariant Faddeev equation which 
traces the binding mechanism of three quarks in a baryon 
to its quark-quark correlations. It is illustrated in Fig. 1, 
where we have already anticipated the gluon-exchange 
interaction that will be motivated in Section II B. The 
Faddeev equation reads explicitly: 

^a.p^s{p,q,P)^ I \K^nLy-^o.p'^'s{p^^\q'^^\P) + 



K. 



(2) 



K 



(3) ^ I Ki) „Vi) p] 

aa'l3P' ^a'P'lSKP ,q 



(3) „(3) 



(7) 



Here, K^'^'' abbreviates the renormalization-group invari- 
ant product of two dressed quark propagators and the 
remaining two-quark kernel /C: 



K 



(a) 



r q(b) q{c) 



(8) 



FIG. 1. (Color online) The covariant Faddeev equation for a baryon amplitude, Eq. (7), in rainbow-ladder truncation. 



where {a, 6, c} is an even permutation of {1,2,3}. The 
covariant nucleon ampHtude ^' carries three spinor in- 
dices {a, /3, 7} for the valence quarks and one index 6 
for the spin- 1/2 nucleon. It depends on two relative mo- 
menta p and q and the total nucleon momentum P, where 
— — is fixed. For details on the kinematics and 
notation, as well as the solution method, we refer the 
reader to Refs. [16, 18, 19]. 

It is noteworthy that Poincare covariance allows for a 
rich spin and quark orbital angular-momentum structure 
in the three-quark Faddeev amplitude. Specifically, the 
nucleon's bound-state amplitude consists of 64 covariant 
basis elements which can be classified in sets of s— , p— 
and d— waves in the nucleon's rest frame [16]. While the 
nucleon ground state is dominated by s— waves, a gluon- 
exchange interaction generates ~ 30% p— wave admixture 
in the nucleon's canonical normalization, and d— waves 
contribute roughly 1% [18]. 

B. Quark propagator and qq kernel 

The kernel (8) of the covariant Faddeev equation in- 
cludes the two crucial ingredients that relate the proper- 
ties of hadrons to the underlying structure of QCD: the 
dressed quark propagator S{p) and the two-quark kernel 
/C which is irreducible with respect to a qq pair. Re- 
peating the steps (4-6) in the quark-antiquark channel 
yields the homogeneous Bethe-Salpeter equation (BSE) 
for a meson that involves the quark-antiquark analogue of 
the kernel IC. The same interaction thereby enters both 
qq and qqq bound-state equations which emphasizes the 
close connection between meson and baryon properties 
in our approach. 

To ensure the pion's nature as the Goldstone boson of 
spontaneous chiral symmetry breaking, the axial-vector 
Ward-Takahashi identity (AXWTI) must be satisfied. 
The AXWTI entails a massless pion in the chiral limit 
and leads to a generalized Gell-Mann-Oakes-Renner re- 
lation [20, 21], and it poses a constraint on the construc- 
tion of the kernel JC by relating it with the kernel of the 
quark DSE, cf. Eq. (40). The simplest kernel to satisfy 
that constraint is the rainbow-ladder truncation which 
amounts to 

/Caa'/3/3' = ^2 p 7aa' 7/3/3', (9) 

where Tj^" = 5^^^ — k^k" /k"^ is a transverse projector 
with respect to the gluon momentum k and Z2 is the 



quark renormalization constant. Eq. (9) describes an it- 
erated dressed-gluon exchange between quark and an- 
tiquark that retains only the vector part ^ 7^ of the 
quark-gluon vertex. Its non-perturbative dressing, to- 
gether with that of the gluon propagator, is absorbed 
into an effective coupling a(fc^) which is modeled. 

There have been various attempts to go beyond the 
rainbow-ladder truncation by employing different strate- 
gies. While Refs. [22-26] analyzed infrared and pion- 
cloud contributions in the quark-gluon interaction di- 
rectly from the vertex DSE, Refs. [11, 27, 28] use a 
Ward-Takahashi identity for the Bethe-Salpeter kernel 
as a guidance for vertex models. These approaches have 
been quite successful in the light meson sector but their 
implementation for baryons is still hampered by techni- 
cal difficulties. For this reason we restrict ourselves to 
the well-established rainbow- ladder kernel of Eq. (9). Its 
main drawback appears in the form of missing pion-cloud 
contributions which are important for the form factors in 
the chiral and low-momentum region [18, 29]. This issue 
is discussed in more detail in section IV. 

The second ingredient in the Faddeev kernel is the 
dressed quark propagator. It is expressed in terms of 
two scalar functions, the quark wave-function renormal- 
ization 1/A{p'^) and the quark mass function M(p^): 

S-\p) = A{p'){z^ + M{p^)). (10) 

The quark propagator satisfies the quark DSE whose in- 
teraction kernel includes the dressed gluon propagator as 
well as one bare and one dressed quark-gluon vertex. In 
rainbow-ladder truncation that kernel becomes identical 
to Eq. (9) and the quark DSE reads 

^apiP) = ^2 {i^ + mo)^p + J ICaa'f3'l3 Sa'l3'{q) ■ (H) 

The bare current-quark mass mo is related to the rcnor- 
malized current mass ruq via 

Z2TO0 = Z2Z„imq ^ ZiTUq , (12) 

where the quark-mass and wave- function renormalization 
constants Z^n and Z2 are determined self-consistently in 
the process of solving the quark DSE. The current-quark 
mass constitutes an input of the equation and can be 
readily varied from the chiral limit up to the heavy-quark 
regime. 

Dynamical chiral symmetry breaking becomes mani- 
fest in the quark DSE if the kernel K. supplies suffi- 
cient interaction strength. The consequence is a non- 
perturbative enhancement of the quark mass function 
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M{p'^) at small momenta that indicates the dynamical 
generation of a constituent-quark mass scale. In prin- 
ciple, such strength would be generated through a self- 
consistent DSE solution for the gluon propagator and 
quark-giuon vertex that enter Eq. (11); see, e.g., [10, 11] 
and references therein. In rainbow-ladder that effect is 
provided by the effective coupling a(fc^) which we choose 
from Ref. [30] as^ 



TTT] 



fc2 
A2 



(13) 



The second term is only relevant at large gluon mo- 
menta where it dominates and is constrained by pertur- 
bative QCD, i.e., it preserves the one-loop renormaliza- 
tion group behavior of QCD for solutions of the quark 
DSE: 



auv(fc^) = 



27r7,.(l-e-^-V^?) 
ln[e2-l + (l + fcVA|cD)'] 



(14) 



with At = 1 GeV, Aqcd = 0.234 GeV, and 7™ = 12/25. 
The first term provides the necessary strength at small 
and intermediate momenta that triggers the transition 
from a current-quark to a dynamically generated con- 
stituent quark. It is characterized by two parameters^: 
an infrared scale A that represents the scale of dynami- 
cal chiral symmetry breaking, and a dimensionless width 
parameter rj. 

In combination with the interaction of Eq. (13), the 
rainbow-ladder truncation has been frequently used in 
Dyson-Schwinger studies of hadrons. By setting the scale 
A via the experimental pion decay constant, the ap- 
proach describes pseudoscalar-meson, vector-meson, nu- 
cleon and A ground-state observables reasonably well, 
see e.g. [32-35] and references therein. In addition, their 
properties are insensitive to a variation of the infrared 
shape of the coupling [30, 34] which is controlled by the 
parameter rj. 



C. Nucleon current 

After having specified the basic ingredients that enter 
the description of baryons in the bound-state approach, 
we turn to the question of resolving the nucleon matrix 
elements of the currents in Eq. (1) in terms of QCD's 
non-perturbative Green functions. A systematic con- 
struction of the nucleon's axial and pseudoscalar current 
in the three-quark framework is provided by the 'gaug- 
ing of equations' method of Refs. [36-38]. It was pre- 
viously applied to derive the electromagnetic current in 



^ We do not expect much difference in our results from a recently 

suggested update [31] of Eq. (13). 
^ The relationship with the parameters {uj, D} used in Ref [30] is 

uiD = A^, r] = A/uj. 



the quark-diquark model [14, 39] as well as in the present 
three-quark setup [18], and a recent computation of the 
pseudoscalar iV — A transition current is based on that 
procedure as well [40]. Moreover, the approach can be 
generalized to yield a quark-level description of general- 
ized parton distributions [41] and a variety of hadronic 
scattering processes [42] . Since the method is not partial 
to the type of the involved current, we can simply adopt 
the derivation in Ref. [18] to the case investigated here. 

The coupling of an external current to a baryon 
amounts to 'gauging' the qqq scattering matrix T. That 
operation acts as a derivative, i.e., it is linear and satisfies 
Leibniz' rule, and we formally denote it by a superscript 
[fj]. The current J 1^1 is the residue of the gauged scat- 
tering matrix T^^l at the hadron's bound-state pole: 



{Pf + M^){P^ + AP) ' 



(15) 



where = '^{pi,qi,Pi) and */ = •^{pf,qf,Pf) are 
in- and outgoing baryon amplitudes with different mo- 
mentum dependencies. From the derivative property 
= -T(T-i)'^'r, in combination with Eqs. (4-5), 
one derives the general expression for a baryon's non- 
perturbative current. It is the gauged inverse T-matrix 
sandwiched between the bound-state amplitudes: 



In the second step we have exploited Eq. (4), 



(16) 



(T-i 



K- 



Go 



= Go (Go 1) Go-K-^ K^^^ K-^ , 



and defined the 'three-body' vertex 

:= (Go^)'^' = (5-1 §5 5-1® 5-1) 



(17) 



(18) 



as the gauged disconnected inverse propagator product. 
To arrive at Eq. (16), we have eliminated the inverse 
kernels that appear in Eq. (17) by inserting the bound- 
state equations K'^^i = Go^i and ^fK'^ = */Go. 
The resulting expression for the current matrix element 
involves the same elements that appear in the bound- 
state equation (6), namely the dressed quark propagator 
and the three-quark kernel K. 

The question remains what 'gauging' of S^^ and K 
means at the quark level. Microscopically, the axialvec- 
tor and pseudoscalar couplings to the quark are repre- 
sented by the respective quark-antiquark vertices pl'^l. 
We can identify them with the gauged inverse quark 

propagator, i.e., (S^^)'''' ~ fI''', and correspondingly: 
51^1 = -S'FI'^1 S. The gauged three-body kernel Xl^l is 
obtained by applying the Leibniz rule to Eq. (6). In 
rainbow-ladder truncation, where the three-body irre- 
ducible contribution to K is neglected and the two-body 
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FIG. 2. (Color online) The two types of diagrams (modulo 
permutations) which contribute to the nucleon's three-body 
current in rainbow-ladder truncation, Eqs. (16) and (19). 



kernel is modeled by gluon exchange, only the direct cou- 
plings to the quarks remain. Then, the ingredients of 
Eq. (16) become 



a=l 
3 



(19) 



where the quark labels {a, 6, c} are again an even per- 
mutation of {1, 2, 3}. The resulting current is illustrated 
in Fig. 2 and consists of an impulse-approximation dia- 
gram and another contribution that involves the gluon- 
exchange kernel. These diagrams are worked out in detail 
in App. B. 

Here we find another practical advantage of the 
rainbow-ladder truncation: due to the structure of 
Eqs. (16) and (19), all ingredients of the nucleon's cur- 
rent are already specified, namely by the ansatz (9) alone. 
This is also true for the dressed axialvector and pseu- 
doscalar vertices that appear in Fig. 2, and we will dis- 
cuss them in the following subsection. 



D. Axial and pseudoscalar vertices 

In order to compute the form factor diagrams of Fig. 2 
we need to specify the dressed pseudoscalar and axi- 
alvector quark-antiquark vertices, i.e., the 'gauged' quark 
propagators, that enter the diagrams. Their proper- 
ties have been established and discussed in detail in 
Refs. [11, 20, 43-45]. Since they have a direct correspon- 
dence with the nucleon's axial and pseudoscalar form fac- 
tors, we find it useful to recollect some of these features 
here and also draw attention to a few properties that are 
perhaps less well-known. 

We can express the coupling of a dressed quark to a 
current with pseudoscalar, vector or axialvector tensor 
structure Fq, with 



Fo e { ^4^75 , Z2ir , Z2l^Y } , 



(20) 



by the quark-antiquark vertex F^^l , where the gauging 
index characterizes the type of Fg . The vertex depends 
on two momenta, the relative qq momentum k and the 



total momentum Q (see Fig. 3). We consider isotriplet 
currents with flavor matrices akl'2., cf. Eq. (1). The ver- 
tex fI^I can be defined as the contraction of the quark- 
antiquark four-point function G with Fq: 



fI^I =GoiGFo = Fo+TGoFo. 



(21) 



where we have exploited Eq. (2) in the quark-antiquark 
case, and Go = S ® S is here the product of two dressed 
quark propagators. This relation is illustrated in Fig. 3. 
The scattering equation (4) immediately yields an inho- 
mogeneous BSE for fI'^1 



fI'^1 = Fo +if GoFt'^1 



(22) 



which allows to compute the vertex for a given kernel K. 
The quark-antiquark scattering matrix T incorporates 



all meson poles, and on the mass shell Q — 
meson with mass ttim it assumes the form 



T 



for a 



(23) 



where 5'm defines the meson's homogeneous bound-state 
amplitude. The central feature that will have an impact 
on the momentum structure of hadron form factors is 
the fact that these poles, via Eq. (21), also appear in the 
vertex fI^I, cf. Fig. 3: 



pip] 



M • 



(24) 



This is true as long as the bound-state wave function 
Go^M has a component in the direction of Fq, which is 
expressed by the residue 



rM[Fo] =Tr / *m Go Fq 



(25) 



For example, the pseudoscalar vertex F5 contains all 
pseudoscalar meson poles, and the quark-photon (i.e., 
vector) vertex F'' inherits the vector-meson poles. The 
respective residues for the pion and the p— meson are 
given by [20, 43] 



Fq = ^4*75 

Fo = Z2ir 



2m„ 



[ro] = - fpmpTl 



(26) 
(27) 



where and fp are the pion and p— meson leptonic de- 
cay constants and niq is the renormalized current-quark 
mass. 

Turning to the axialvector vertex F5, it is advanta- 
geous to write the vertex as a sum of purely transverse 
and longitudinal contributions with respect to the total 
momentum Q: 



F^(fc, Q) = F^^(fc, Q) + ^ F5.L(fc, Q) , (28) 




■-- + 




-mil 




FIG. 3. (Color online) Diagrammatical representation of Eq. (21) and its onshell limit, Eqs. (23-25). 



SO that Fg — F5 Analyticity at = implies cor- 
relations between the transverse and longitudinal parts, 
see Eq. (37) below. Since axialvector mesons are trans- 
verse, their bound-state poles will appear only in Fg^ 
whereas the longitudinal part T^ j^ contains pseudoscalar 
poles. This can be seen from inserting 



757 



757t 



■75S 



(29) 



together with the decomposition (28) into Eq. (21), upon 
which the equation decouples: 



(l + rGo)(Z2757T) 

r5,L = (i + rGo)(Z275(^). 



(30) 



From Eq. (25) one infers that the longitudinal tensor 
component 75^ can overlap with pseudoscalar pole struc- 
tures encoded in T, for example the pion's bound-state 
wave function, and that overlap is just the definition of 
the pion decay constant [20]: 



To = ^275^ 



r^r [Eo] 



(31) 



Thus, Fs^i will exhibit pseudoscalar poles at timelike val- 
ues of the total momentum-squared , and axialvector 
poles will appear in Fg^. Similarly, the inhomogeneous 
axialvector BSEs (22) for those two vertices decouple as 
well, and the longitudinal equation is identical to the 
pseudoscalar BSE except for the different driving term 
■^275^ instead of .^4175. The pole behavior will there- 
fore be recovered in the solution of the inhomogeneous 
BSEs, and by implementing these vertices in the axial 
current of the nucleon, it will translate to the nucleon's 
axial form factors in the timelike region, cf. Section IIF 
The precise relation between the pseudoscalar vertex 
and the longitudinal part of the axialvector vertex is ex- 
pressed by the AXWTI which reads in the flavor-triplet 
case: 



r5,L(fc,g) + 2m,F5(fc,Q) = 



(32) 



where k± — k±Q/2 are the two quark momenta. Eq. (32) 
entails that the pseudoscalar poles on the left-hand side 
must compensate each other, and that the longitudinal 
part of the axialvector vertex is completely specified by 
the pseudoscalar vertex and the dressed quark propaga- 
tor alone. 



The AXWTI can be made explicit in a given basis 
decomposition of the vertices. The transverse part of the 
axialvector vertex involves eight tensor structures which 
can be chosen as [46]: 



= 75 



7T(/i+*/2fc-0^) + 

+ ki^k-Q{ih + + 
+ k^ {f7^ + if8k-Qh(^) 



(33) 



where the fi{k'^, k ■ Q, Q^) are scalar dressing functions 
and 7^, fc^ are transverse with respect to the photon mo- 
mentum Q. The angular prefactors k-Q were attached to 
guarantee positive charge-conjugation parity (for quarks 
with equal mass), so that the vertex carries the quantum 
numbers J^^ = 1++. The factors i ensure that all dress- 
ing functions are real if fc^ G R+ and G R. Likewise, 
the longitudinal part can be decomposed as 



F 



5,L 



«75 



(34) 



and an analogous decomposition holds for the pseu- 
doscalar ^ vertex: 



= «75 



(35) 



The AXWTI (32) yields the following relations between 
the components of T^j^ and F5 which are valid for all 
current-quark masses: 



gi + 2mq hi =2T.B , 
g2 + 2mq /i2 = -Sa , 
5(3 -I- 2mq /13 = -2 , 
34 -I- 2mq /14 = . 



(36) 



Here, A{k'^) and B{k'^) = M{P)A{P) are the dress- 
ing functions of the inverse quark propagator S^^{k) = 
ifi: A{k'^) + B{k^), and we used the abbreviations 



:= 



F{kX)+F{kl) 



Ap := 



F{k\)-F{kl) 
kl-kl 



with F e {A,B}. 

In the limit — > 0, regularity leads to an addi- 
tional relation between the longitudinal and the trans- 
verse dressing functions in Eqs. (33-34): 



51 = 52 + /i = 53 + /? = 34 + /a = , 



(37) 
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and thereby correlates four of the eight dressing functions 
of the axial- vector vertex with those of the pseudoscalar 
vertex via Eq. (36). For instance, one obtains at — > 0: 
/ii(fc2) = B{P)/mq. With Eqs. (24) and (26), the chiral- 
limit behavior of the pseudoscalar vertex then becomes 



m„ 



2m„ 



(38) 



where ^'tt is the pion's canonically normalized bound- 
state amplitude, with the same structure as Eq. (35) ex- 
cept for the flavor convention cr^ instead of (7^/2. This 
yields the well-known chiral-limit relation between the 
dominant pion dressing function and the scalar dressing 
of the quark propagator [20] : 



/.W(fc2) 



(39) 



These considerations finally provide a quick way to 
check the Gell-Mann-Oakes-Renner relation which is al- 
ready implicit in Eq. (26). If Eq. (31) defines the pion 
decay constant, equating the residues of the vertices T^^l 
and Ts in the AXWTI (32) yields the quoted value for 

in Eq. (26). On the other hand, the trace with ^4^75 
in (26) filters out the dominant component of the pion's 
wave function Go^'tt which, in the chiral limit, satisfies 
a relation analogous to Eq. (39), except that the inverse 
dressing B{k^) is replaced by the scalar dressing of the 
quark propagator itself. Its integrated version is just 
the quark condensate, i.e., the trace over the chiral-limit 
quark propagator. 

While the properties discussed so far are exact in QCD, 
a practical calculation of the vertices and F5 from the 
inhomogeneous BSE (22) requires an ansatz for the kernel 
K to operate with. Since the characteristics of the pion 
as the Goldstone boson of spontaneous chiral symmetry 
breaking are intimately related with the AXWTI, it is 
vital that any kernel ansatz complies with that relation. 
Through the inhomogeneous BSE (22), the AXWTI can 
be reformulated as a relation between the dressed quark 
propagator and the qq kernel K: 

k' (40) 
= - «75 + «75 S(fc_)]„^ , 

where 5](fc) is the quark self-energy (not to be confused 
with the function from above) that appears in the 
quark DSE (11): S-\k) = Z2{Ht + mo) + ^{k). Eq. (40) 
and analogues thereof can be used as a construction prin- 
ciple for suitable kernel ansatze [47, 48], and it is not 
difhcult to prove that the relation is preserved by the 
rainbow-ladder kernel of Eq. (9). 

We further recall that similar considerations apply for 
the quark-photon vertex F'' whose longitudinal part, de- 
fined in analogy to Eq. (28), must reproduce the electro- 
magnetic Ward-Takahashi identity: 

Q^T^ik, Q) = TL{k, Q) = S-\k+) - S-\k^) . (41) 



Thus, the longitudinal part of the vector vertex is com- 
pletely fixed by electromagnetic gauge invariance, i.e. by 
the dressed quark propagator alone. Combined with the 
analyticity requirement at = 0, Eq. (41) leads to the 
Ball-Chiu construction for the vertex [46] , augmented by 
a purely transverse term that includes an array of vector- 
meson poles. Upon implementation in the form factor 
diagrams, the latter is responsible for the timelike pole 
structure in electromagnetic form factors and the under- 
lying reason for vector-meson dominance. Similarly to 
Eq. (40), the electromagnetic WTI can be formulated 
as a relation between the qq kernel and the quark self- 
energy: 

J^K^^,p,4Sik'_,)-Sik'_)]^,^, ^^^^ 

-[I](fc+)-S(fc_)]„^ , 

which is also satisfied by the rainbow-ladder kernel and 
consequently by the vector vertex as a solution of its 
inhomogeneous BSE. 

We have now collected all ingredients for solving the in- 
homogeneous Bethe-Salpeter equations for the axialvec- 
tor and pseudoscalar vertices. They take the explicit 
form 

F^_^(fc,Q) = Z2 757^-^KT], 

r5.L(fc,Q) = ^2 75 ^-2:[F5^L] , (43) 

^5{k, Q) = Zi 175 - 1 [F5] , 
where the integrals on the right-hand side are given by 



'lap [F] 



[S{k'_,)T{k',Q)S{k'_)\^ , , (44) 



with the rainbow-ladder kernel K, from Eq. (9) . Since the 
rainbow-ladder truncation obeys the AXWTI per con- 
struction, the relations (36-37) for the resulting vertex 
dressing functions are numerically satisfied to high accu- 
racy. With the axialvector and pseudoscalar vertices ob- 
tained from Eqs. (43), all building blocks of the nucleon 
current matrix element in Eq. (16) are now specified. 

To summarize, the gauging method provides in com- 
bination with Eq. (21) a systematic and intuitive way of 
computing hadronic matrix elements of local quark cur- 
rents. If a current Fq probes a hadron, it does so via 
Eq. (16), i.e., the interaction is resolved in a coupling 
to the quarks and the kernels (which are microscopically 
again resolved into a coupling to quarks). The interac- 
tion of the current with the dressed quark is described 
by Eq. (21) which amounts to the sum of a pointlike part 
plus all possible reaction mechanisms between quark and 
antiquark which constitute the T-matrix. In practice, 
the latter part is accessible by resumming the full ver- 
tex via solving its inhomogeneous BSE. We finally note 
that Eq. (16) is valid for arbitrary kernel ansatze. Thus, 
if the two- and three-body kernels are known diagram- 
matically, the method can be applied for implementing 
interactions beyond rainbow-ladder as well. 
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III. AXIAL AND PSEUDOSCALAR FORM 
FACTORS 

We continue by discussing the implications of the re- 
lations collected in Section II D for the axialvector and 
pseudoscalar current matrix elements in the nucleon. 
The axialvector current is characterized by two form 
factors: the axial form factor Ga{Q^) whose value at 
zero momentum transfer is the nucleon's axial coupling 
constant and the induced pseudoscalar form factor 
Gp{Q^). They constitute the nucleon matrix elements 
of the axialvector current in Eq. (1) via 



J^{P,Q)^^ilAGAr + G 



2M 



(45) 



where Pi and P/ are the incoming and outgoing nu- 
cleon momenta, Q ^ Pf — Pi is the momentum transfer, 
P = {Pi + Pf)/2 is the average momentum, and M is the 
nucleon mass. A further possible structure proportional 
to P'^ is forbidden by charge-conjugation invariance. In- 
stead of nucleon spinors we use the positive-energy pro- 
jectors A:^'' = {1 + |'/,i)/2, so that Eq. (45) is a matrix 
in spinor space. Similarly, the pseudoscalar current is 
parametrized by the pseudoscalar form factor Gz{Q'^): 



MP.,Q)^G5Ai 175 A^^ 



(46) 



It is instructive to split the axialvector current (45) into 
purely transverse and longitudinal contributions with re- 
spect to the momentum transfer Q. Inserting the trans- 
verse 7— matrix from Eq. (29), together with the relation 
A^75(^A; = -2iMA^75 AV, yields: 



J^{P,Q) ^Aij5 {Ga1^~Gl 



2Mt 



A^ 



+ ' 



where t — Q"^ / {ANP) . Gp and Gl are related via 



Gp — - {Ga — Gl) 

T 



(47) 



(48) 



and analyticity for Gp ai vanishing momentum transfer 
requires Gl{0) = Ga{0). 

The relations of Section II D directly translate to the 
nucleon level since the form-factor diagrams of Fig. 2 are 
identical in the axialvector and pseudoscalar cases except 
for the type of vertex that is involved. Thus, the proper- 
ties of the qq vertices T'^j,, T^^l and Tr, are reflected in 
the form factors Ga, Gl and G5, respectively. For time- 
like values < 0, the axial form factor Ga must inherit 
the transverse axial- vector poles encompassed in F^^, 
whereas the pseudoscalar pole structure will appear in 
both Gl and G5 . Eq. (48) implies that the induced pseu- 
doscalar form factor Gp will inherit both pseudoscalar 
and axial- vector poles. 

On the pion's mass shell, the residue of the pseu- 
doscalar form factor G5 is the pion-nucleon coupling con- 
stant which can be made explicit by defining 



G5(g2 



Q2 + fjil niq 



GttNN 



(49) 



with G^Ar7v(— TO^) = g-nNN- This is the analogue of the 
vertex relations (24) and (26) at the form-factor level, 
and the missing factor 2 is a consequence of the flavor 
conventions: the qq vertices carry flavor matrices crfc/2 
and the pion's bound-state amplitude a factor ak- 
The AXWTI (32) at the nucleon level becomes 



g^Jg^ + 2mg J5 = 



(50) 



and thereby relates the longitudinal with the pseu- 
doscalar form factor [15, 49]: 



^ M ^ g2 + ml M 



Gjr 



NN ■ 



(51) 



Combined with the analyticity requirement for Gp from 
Eq. (48), this yields the Goldberger-Treiman relation: 



Ga(0) 



M 



GttNn{^) 



(52) 



which is valid for all current-quark masses. Thus, the 
axial and pseudoscalar currents (45-46) are described by 
two independent form factors Ga and Gt^nn which, at 
— 0, are related via the Goldberger-Treiman relation. 
Implementing the AXWTI and the analyticity con- 
straint explicitly, one can further define 



Ga{Q^) = ^ G7rAfjv(g^) 



4Af2 



P(Q') 



(53) 



so that the quantity R{Q^) represents the deviation from 
the Goldberger-Treiman relation at non-zero momentum 
transfer. In that way the two currents are parametrized 
by two independent form factors Ga{Q'^) and R{Q^), or 
G-nNNiQ'^) and R{Q^). This yields for the induced pseu- 
doscalar form factor via Eqs. (48) and (51): 



Gi 



AJVPGa -mlR AMU 



NN 



R. 



(54) 



If Gp is expressed through Ga, the form factor R enters 
in combination with ml; if Gp is related to G-^nn, it 
appears as the non-resonant remainder at the pion pole. 
The pion-pole dominance assumption which corresponds 
to i? = in Eq. (54) is therefore justified, at least in the 
low quark-mass region and the neighborhood of the pion 
pole. The contribution from R is however not necessarily 
small in Eq. (53) itself since that would entail the validity 
of the Goldberger-Treiman relation for all Q^. 

We finally recall that the 'Goldberger-Treiman discrep- 
ancy' in the usual sense is defined as 



Agt = 1 - 



Ga{0) 

JT Gr,NNi~ml) 



(■52) ^ _ 



G. 



NN 



^ 

G-^NNi-ml) 



(55) 



and measures the discrepancy of G^^nn between — 
and the pion pole at a given current-quark mass, and 
thereby the distance from the chiral limit. 
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FIG. 4. (Color online) Results for the nucleon's axial charge qa and axial radius ta as a function of the pion mass. We compare 
to lattice results from Refs. [6-8, 50] and the chiral expansion of Ref. [5] that includes explicit A degrees of freedom. Stars 
denote the experimental values. 



IV. RESULTS AND DISCUSSION 



Before discussing the results, let us briefly revisit the 
successive steps that are performed in the present calcu- 
lation. Once the rainbow-ladder kernel is specified via 
Eqs. (9) and (13), we solve the Dyson-Schwinger equa- 
tion (11) for the quark propagator, the covariant Fad- 
deev equation (7) for the nucleon amplitude, and the in- 
homogeneous Bethe-Salpeter equations (43-44) for the 
axialvector and pseudoscalar quark-antiquark vertices. 
These quantities enter the nucleon's current matrix el- 
ements (16) whose decomposition is illustrated in Fig. 2. 
The explicit construction of the form-factor integrals is 
discussed in detail in App. B. The axial, longitudinal and 
pseudoscalar form factors Ga{Q^), Gl{Q^) and G^{Q'^) 
are finally extracted from the Dirac traces in Eq. (B16). 

Within the rainbow-ladder truncation, the only model 
input in our calculation is the effective quark-gluon in- 
teraction from Eq. (13). It depends on two infrared pa- 
rameters, the scale A and the width parameter rj. The 
calculation is carried out with fixed A = 0.74 GeV that 
is chosen to reproduce the experimental pion decay con- 
stant at the physical pion mass. The width parameter 
rj = 1.8±0.2 is a measure of the model uncertainty and 
leads to the colored bands in Figs. (4-6). By modifying 
the current-quark mass that enters the quark DSE as an 
input, we can further study the quark-mass dependence 
of the form factors. The quark mass can be mapped onto 
the pion mass by simultaneously solving the BSE for the 
pion with the same input. The resulting values for the 
pion decay constant, the nucleon mass, and various other 
observables related to the axial and pseudoscalar currents 
are tabulated in Table I for three different pion masses. 



A. Quark-mass dependence of qa 

Our result for the axial charge gA = G^(0) as a func- 
tion of the squared pion mass is shown in the left panel 
of Fig. 4. Its experimental value gA = 1.2695(29) is pre- 
cisely known from neutron /3 decay [51]. Depending on 
the width parameter in the quark-gluon interaction, 
our calculated result gA — 0.99(2) at the physical u/d- 
quark mass underestimates that value by 20 — 25%. It 
also falls below recent lattice data but slowly increases 
with the pion mass and approaches the chiral expansion 
from Ref. [5] at higher quark masses. 

The determination of gA in chiral perturbation the- 
ory has been notoriously difficult due to strong cancela- 
tions at leading and next-to-leading order in the chiral 
expansion [1, 4]. An important role in the nucleon's ax- 
ial structure has been attributed to the A (1232) [52-54]. 
The implementation of pion loops with internal A de- 
grees of freedom in the chiral expansion produces indeed 
a relatively fiat pion-mass dependence of gA [5] which is 
plotted in Fig. 4 for comparison. Lattice results for gA 
with pion masses below 300 MeV are now available [6-8] 
but have not been able to clarify the issue as they still 
underestimate the axial charge by 10 — 15%, cf. Fig. 4. 
They are, however, strongly sensitive to finite- volume ef- 
fects which has been interpreted as a signal of missing 
pion-cloud contributions on the lattice [55] . 

A similar interpretation arises in the Dyson-Schwinger 
approach. From a microscopic point of view, pion ex- 
change corresponds to a gluon resummation. While the 
pion's bound-state pole is recovered in the qq scatter- 
ing matrix in rainbow-ladder truncation, which corre- 
sponds to a gluon resummation in the s-channel, pion- 
cloud contributions are generated by t— channel pion ex- 
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Mn 


gA 




mA/MN 




GttNn{Q) 


g-nNN 


Agt 


Exp. 


0.138 


0.092 


0.94 


1.27 


1.0. ..1.3 


1.1... 1.4 


1.32(3) 


12.9 


13.2 


0.02 




0.138 


0.092 


0.94(1) 


0.99(2) 


1.28(6) 


1.36(7) 


0.95 


10.2(3) 


10.4(3) 


0.01 


Calc. 


0.522 


0.115 


1.25(2) 


1.05(1) 


1.51(4) 


1.21(5) 


0.89 


11.5(2) 


13.4(5) 


0.14(1) 




0.761 


0.134 


1.51(3) 


1.10 


1.69(3) 


1.12(4) 


0.86 


12.6(2) 


16.4(6) 


0.23(2) 



TABLE I. Results for various observables, as discussed in the text, at three different pion masses and compared to experiment. 
The input scale is chosen such that the pion decay constant reproduces the experimental value at = 138 MeV. The 
parentheses indicate the dependence on the infrared parameter rj of Eq. (13). m-^, Mn and the axial mass tua are given in 
GeV, the remaining quantities are dimensionless. 



change [26] whose relevant gluon topologies are not cap- 
tured in rainbow-ladder. In fact, missing pion-cloud 
contributions have turned out to be the main source 
of discrepancy in the chiral and low-momentum struc- 
ture of the nucleon's electromagnetic form factors com- 
puted in the present setup [18]. For example, the nu- 
cleon's isoscalar anomalous magnetic moment Ks, where 
leading-order chiral corrections cancel, is accurately re- 
produced by the Faddeev calculation; the nucleon charge 
radii underestimate the experimental values but converge 
with lattice data at larger quark masses; and the low- 
behavior of the form factors shows missing structure 
whereas one finds reasonable agreement with experiment 
at larger Q^. These observations suggest to identify the 
rainbow-ladder truncated nucleon with the 'quark core' 
in chiral effective field theories. In the case of axial form 
factors it would imply that the discrepancy between our 
result for gA and its experimental value (and its chiral 
expansion) mainly owes to chiral cloud effects as well. 

Irrespective of that, we find that the Goldberger- 
Treiman relation of Eq. (52) is satisfied quite accurately. 
The relation 

is realized at the permille level at the physical u/d mass 
and still reproduced up to 1 — 2% at the largest quark 
mass where our numerical treatment becomes increas- 
ingly inaccurate. We note that a translationally invari- 
ant cutoff regularization in the quark DSE and vertex 
equations [44] turned out to be crucial in order to en- 
sure its validity; a hard numerical integral cutoff leads 
to violations of several percent in Eqs. (37) and (56). 
This procedure also required a slight readjustment of the 
input scale A compared to Ref. [18] to reproduce the ex- 
perimental pion decay constant. 

As we have discussed in Sections II D and III, Eq. (56) 
is a direct consequence of the underlying dynamics. The 
rainbow-ladder kernel satisfies the AXWTI (32) by con- 
struction, and that property translates to the quark 
propagator, the quark-antiquark vertices, the nucleon 
bound-state amplitude and finally also to the form fac- 
tor level. The Goldberger-Treiman relation entails that 
the same underestimation that is visible in Fig. 4 for 
the axial charge also occurs in G7rArAr(0). The experi- 
mental value for the pion-nucleon coupling constant is 



g-KNN — 13-2 [56], and from the Goldberger-Treiman re- 
lation one obtains the value G'^ArAr(O) = 12.9 at vanish- 
ing momentum transfer. Our result GT^NNiO) = 10.2(3) 
underestimates that phenomenological value again by 
20 - 25%, cf. Table I. 



B. Momentum dependence of Ga 

While the axial coupling constant is determined from (3 
decay, experimental data for the momentum dependence 
of Ga{Q^) have been obtained in quasielastic neutrino 
scattering off nucleons or nuclei and charged pion elec- 
troproduction [1]. The existing data can be parametrized 
by a dipole ansatz, 



Ga(Q') 



9A 



{l + Qym\Y 



(57) 



which determines the axial mass mA- The axial radius 
ta is obtained from the slope of the axial form factor at 
g2-0: 



GAiQ') 



= 1 - 



= -6 



G^AiO) 



9A 



9A 



(58) 



Inserting the dipole form yields ra\ — 12 /r\, which can 
also be used to define the axial mass independently of 
the dipole parametrization. 

Due to model-dependent extractions of Ga from the 
respective cross sections, the experimental value of the 
axial mass is not well constrained. Pion electroproduc- 
tion and older neutrino scattering experiments yield val- 
ues around niA ^ 1 GeV [1, 58] whereas recent data from 
MiniBooNE and K2K favor higher central values up to 
niA ^ 1.3 GeV [59, 60]. The origin of this discrepancy is 
unclear and could be a consequence of nuclear medium 
effects [61] or a deviation from the dipole form [62]. 

In the right panel of Fig. 4 we show our result for 
the quark-mass dependence of the axial radius ta and 
compare it to lattice data. The corresponding results for 
the axial mass niA = Vl2/rA are listed in Table I for 
three different pion masses. In contrast to gA, we find 
reasonable agreement with the lattice data throughout 
the current-mass range. The results are also compatible 
with experiment at the physical u/d mass as long as a 
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FIG. 5. (Color online) Results for Ga{Q^) and G',rJVjv(Q^) at the physical u/d mass. We compare to a selection of lattice 
results at different pion masses [7, 8, 57]. The dashed line with the hatched area in the plot for Ga is the experimental dipole 
form (57) with an axial mass ttia = 1.15 ± 0.15 GeV. The phenomenological curve for GttNn{Q^) is obtained from the dipole 
ansatz for Ga by assuming the validity of the Goldberger-Treiman relation for all , i.e., Eq. (53) with R — Q. 



somewhat larger experimental value for the axial mass is 
assumed, i.e., uia ^ 1-2 GeV. 

The Q^— evolution of Ga is depicted in Fig. 5 and com- 
pared to lattice results and the phenomenological dipole 
parametrization of Eq. (57). At low Q^, our result under- 
estimates the lattice data and the dipole curve whereas 
they converge above ~ 1 . . . 2 GeV^ . This is again a 
similar behavior as in the case of electromagnetic form 
factors, and it might indicate missing chiral cloud effects 
in the low-momentum region. From that point of view, 
our calculation favors again a somewhat larger experi- 
mental value for the axial dipole mass since for smaller 
values of tua the dipole curve will fall below our result 
already in the low-Q^ region. 

In general, the low-Q^ behavior of the form factors is 
dominated by their timelike singularity structure. This 
feature is quite transparent in the Dyson-Schwinger ap- 
proach where it is induced by meson poles in the vertices 
at timelike which have their origin in the qq scattering 
matrix, cf. Eqs. (23-24). While the relevant mass scale in 
the case of electromagnetic form factors is the p— meson, 
Ga corresponds to the transverse part of the axialvector 
qq vertex F^ ^ which has meson poles in the axialvector 

(J^*^ — 1"*""*") isovector channel. Thus, the axial mass 
will be dominated by the properties of the axialvector 
meson ai(1260) and, to some extent, also by its radial 
excitations. 

While wc cannot directly access the form factors for 
timelike values of without a complex continuation 
in the Faddeev amplitude, computing the qq vertices at 
timclike momenta is unproblematic. From that we can 
conclude that our result for Ga must exhibit a pole at 
= —rn1^_^ . Thus, if Ga were to follow a perfect dipole 



form, our calculated axial mass would coincide with the 
result for obtained in the same setup. The mass of 
the ai meson is considerably underestimated in rainbow- 
ladder truncation, with TOqi — 0.90(1) GeV depending 
on the parameter r?; see Refs. [45, 63] for a comprehen- 
sive collection of ground-state and excited meson masses 
in rainbow-ladder. This deficiency is well known and cor- 
rections beyond rainbow-ladder can be held accountable 
for it [11, 25]. On the other hand, our result for the axial 
mass extracted from the axial radius (58) is = 1.28(6) 
GeV which implies that Ga shows a sizeable deviation 
from the dipole parametrization. 

It is instructive to revisit the current-mass dependence 
of GaIQ"^) in view of these considerations. We note that 
the calculated ratio vria^ /Mjv decreases only slowly with 
the quark mass, cf. Table I. A weak current-mass de- 
pendence of dimensionless ratios is realized for various 
observables in the Dyson-Schwinger approach [64] and 
can be traced back to the interplay of the two scales, 
A and mg, that enter the calculation of hadron proper- 
ties. A represents the scale of dynamical chiral symmetry 
breaking and enters the definition of the rainbow-ladder 
interaction (13). It controls the chiral regime and, with 
the exception of the pion mass, entails a chiral-limit scal- 
ing of hadron masses with A [65] . When increasing the 
current-quark mass m^, the impact of dynamical chiral 
symmetry breaking diminishes since it is roughly quark- 
mass independent, and hadron masses scale with vriq. 

A weak dependence of TOq^ /Afjv on the quark mass im- 
plies that the pole position of the ai— meson in the vari- 
able Q^/M^, and therefore the overall shape of the form 
factor Ga as a function of Q"^ /M'^ , will be insensitive to 
a change in the current-quark mass. We find indeed that 
the Q^/M^— behavior of the form factor Ga in Fig. 5 
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FIG. 6. (Color online) Result for the induced pseudoscalar 
form factor Gp{Q^) at the physical u/d mass. The lat- 
tice references are the same as in Fig. 5 for G'a(Q^). The 
phenomenological band corresponds to the pion-pole domi- 
nance assumption, i.e., Eq. (54) with i? = 0, in combination 
with the dipole parametrization for Ga with an axial mass 
ruA = 1.15 ±0.15 GeV. 



does not appreciably change when going to higher quark 
masses. The same can be said about the lattice data 
of Refs. [6, 8]: when plotted over /AP, the spread be- 
tween the data obtained at various pion masses (of which 
only a subset at the lightest available pion masses are 
shown in Fig. 5) is significantly reduced, and they fall at 
a relatively narrow band that agrees well with the dipole 
parametrization at the physical point. Similar properties 
have been observed for various other form factors, e.g., 
the electromagnetic form factors of the nucleon [18] and 
the A resonance [66]. Of course, genuine chiral features 
such as non-analyticities stemming from the pion cloud 
and the opening of decay channels would lead to a devi- 
ation from this behavior; however, such features are not 
captured by our present truncation. 



C. Pseudoscalar form factor Gttnn 

We now turn to the strong pseudoscalar form factor 
G-nNNiQ'^) which is extracted from Gz{Q'^) upon remov- 
ing the pion pole contribution, cf. Eq. (49). The form 
factor is not directly observable except at the onshell 
point = —m'^, where it reproduces the pion-nucleon 
coupling constant g-^NN- 

The observations from the previous subsection can be 
carried over here. The pole structure in the pseu- 
doscalar vertex that enters the form factor diagrams en- 
tails that GttNN is dominated by the properties of the 
first radial excitation of the pion, i.e., the 7r(1300). Its 
mass in rainbow-ladder truncation is 1.1 GeV, with a 



current-mass dependence similar to that of the ai me- 
son [67, 68]. As a consequence, the form factor G-nNN 
as a function of Q"^ /M'^ shows again a similar behavior 
throughout the current-mass range. 

In the right panel Fig. 5 we show the momentum 
dependence of G^^nn at the physical point. The phe- 
nomenological band denotes the continuation of the 
Goldberger-Treiman relation (52) to non-zero using 
the dipole parametrization for Ga- The previous con- 
siderations make clear that there is no reason for the 
Goldberger-Treiman relation, which relates G^^Nj^f and 
Ga at — due to analyticity, to hold at non-zero 
momentum transfer. The properties of both form factors 
are dominated by different (pseudoscalar and axialvec- 
tor) meson poles with no inherent connection. Neverthe- 
less, one would still expect a qualitatively similar space- 
like behavior of the form factors since the masses of the 
ai(1260) and 7r(1300) are sufficiently close and share a 
similar quark-mass dependence. 

Fig. 5 illustrates that our result for the pseudoscalar 
form factor deviates indeed from the Goldberger-Treiman 
relation at > whereas it is compatible with the 
lattice data of Ref. [57]. In Table I we also show the 
result for the pion-nucleon coupling constant g^NN at 
three different quark masses. We obtain it by extrap- 
olating GjrNN from spacelike momenta to — —rn^. 
Since the pion mass vanishes in the chiral limit and fol- 
lows the Gell-Mann-Oakes-Renner relation (mj ~ niq) at 
small quark masses, the location of = will shift 
to larger timelike values upon increasing the quark mass. 
The form factor GT^NNiQ'^ /M^) depends only weakly on 
the quark mass, and therefore g-nNN will grow faster with 
increasing quark mass than the static value G7rArAf(0); see 
Ref. [40] for an analogous discussion of the ANtt coupling 
constant. As a consequence, the Goldberger-Treiman dis- 
crepancy from Eq. (55), which measures the relative dif- 
ference in Gttnn between = and = —itl^, rises 
with the squared pion mass, see Table I. 

D. Induced pseudoscalar form factor Gp 

We conclude our analysis with the induced pseu- 
doscalar form factor Gp{Q^) that enters the axial current 
of Eq. (45). Experimental data for Gp are sparse: its 
Q^— dependence can be extracted from pion electropro- 
duction, whereas ordinary muon capture on the proton 
{fi~ + p ^ Vfj^ + n) determines the induced pseudoscalar 
coupling constant gp: 

5p = ^Gp(Q2 = 0.88M^), (59) 

where is the muon mass. The recent value reported 
by the MuCap Collaboration, gp = 7.3 ± 1.1 [69], is con- 
sistent with chiral perturbation theory but smaller than 
the previous world average [1, 2]. 

From Eqs. (48) and (54) it is clear that Gp{Q'^) is 
dominated by the pion pole at = —rn^. This implies 
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that the pion-pole dominance assumption 



4M^ Ga 



mi. 



(60) 



describes the spaceUke properties of the induced pseu- 
doscalar form factor reasonably well. Such a behavior is 
clearly visible in our result shown in Fig. 6. In partic- 
ular, our underestimation of the axial charge trans- 
lates to the induced pseudoscalar coupling constant and 
yields gp = 6.7(1). Since Gp is the linear combination of 
the axial and longitudinal form factors Ga and Gl via 
Eq. (48), we can also infer information about its struc- 
ture beyond the pion pole: namely, it must include all 
further pseudoscalar (0 ^) and axialvector (l^"*") pole 
structures as well. 

In principle, the longitudinal form factor Gl allows to 
check the vaUdity of the AXWTI (51) that relates Gl 
with the pion-nucleon form factor Gt^nn- Unfortunately 
we find that Gl is rather sensitive to the angular de- 
pendencies in the nucleon's Faddeev amplitude. In the 
present numerical treatment, the nucleon amplitude is 
obtained from solving its Faddeev equation in the rest 
frame. In the form factor calculation, a 'Lorentz boost' 
to a moving frame amounts to a complex continuation 
of the amplitude's dressing functions in its angular vari- 
ables which is usually implemented via Chebyshev ex- 
pansions. Potential convergence problems with growing 

can be alleviated by imposing a Q^-dependent cutoff 
in the Chebyshev expansion [18]. Even then, however, we 
find that the axial and pseudoscalar form factors at larger 

are more sensitive to the numerics than their electro- 
magnetic counterparts. This is especially true, even at 
low momentum transfer, for the longitudinal form fac- 
tor Gl{Q^) which is numerically not well under control. 
While that numerical instability does not affect the be- 
havior of the induced pseudoscalar form factor which is 
dominated by Ga and the pion pole, it impedes a test of 
the AXWTI. 

We want to emphasize that the problem could be en- 
tirely eliminated by solving the Faddeev equation in each 
boosted frame (i.e., at each Q^) anew. That would al- 
low for a clean determination of all three form factors 
Ga, Gp and G^^nn up to ~ 7... 9 GeV^ without 
the need for analytic continuation in the angular vari- 
ables. The upper limit is due to the singularities in 
the quark propagators that appear in the integrands; it 
could be circumvented as well via residue calculus. How- 
ever, these strategies are both numerically and concep- 
tually challenging and beyond the scope of the present 
work. 



as well as their microscopic building blocks. All ingre- 
dients were computed selfconsistently within a rainbow- 
ladder truncation which describes the binding of the nu- 
cleon through iterated dressed gluon exchange between 
the quarks. Since pion-cloud effects are not implemented, 
our results represent the nucleon's quark core. 

We analyzed the properties of the axial (Ga), in- 
duced pseudoscalar (Gp) and pion-nucleon form factor 
{GttNn) from the underlying quark-antiquark vertices 
on which they depend. The spacelike structure of the 
form factors is connected to the properties of axialvec- 
tor and pseudoscalar-meson poles in the timelike region. 
Ga is dominated by the axialvector meson ai(1260) and 
its radial excitations; Gp is governed by the pion pole, 
and Gt^nn is related to the 7r(1300). The Goldberger- 
Treiman relation follows as a consequence of the axialvec- 
tor Ward-Takahashi identity and analyticity which are 
both satisfied at the quark-gluon level. There is however 
no obvious inherent connection between the form factors 
Ga and G-^mn at non-zero momentum transfer. 

The form factor results at = underestimate the 
experimental values by 20 — 25% whereas they are consis- 
tent with the phenomenological dipole form for Ga and 
the pion-pole dominance ansatz for Gp at larger . The 
nucleon's axial charge gA falls below the lattice data in 
the low quark-mass region and approaches the result from 
chiral effective field theory at larger pion masses. These 
features might be signals of missing pion-cloud effects in 
the chiral and low-momentum regions. 

We encountered difficulties in reproducing the axi- 
alvector Ward-Takahashi identity beyond small as it 
turned out to be quite sensitive to the numerical details 
of the Faddeev amplitude. A clean determination of the 
form factors at larger would require a moving-frame 
solution of the Faddeev equation. This is numerically 
challenging and remains a task for future investigations. 
Moreover, in order to decide whether the missing struc- 
ture in the form factors is indeed attributable to pion- 
cloud effects, a consistent description of baryons beyond 
the rainbow-ladder truncation is necessary. Steps in that 
direction are planned. 
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Appendix A: Euclidean conventions 

We work in Euclidean momentum space with the fol- 
lowing conventions: 



p ■ q 



= ^Pkqk, p'^^p-p, ^^p-j. 



(Al) 



fc=i 



A vector p is spacelike if > and timelike if p^ < 0. 
The hermitian 7— matrices 7'' = (7^)^ satisfy the anti- 
commutation relations {7'', 7'^} = 2 6'^", and we define 



7 



(A2) 



In the standard representation one has: 















/ 


,7^ = 


I 






7= . n '7=K. ^ '7' 



where are the three Pauli matrices. The charge con- 
jugation matrix is given by 



C 



4 2 
7 7 , 



c- 



-c, 



(A3) 



and the charge conjugates for (pseudo-) scalar, (axial-) 
vector and tensor amplitudes are defined as 

T{p,p) ■.^cr{-p,-pfc^ , 
r°'{p,P) := -cr°'{-p,-pf , (A4) 
TI^'^{p,P) := Cr"'5(-p,-P)^C^, 

where T denotes a Dirac transpose. Four-momenta are 
conveniently expressed through hyperspherical coordi- 
nates: 



p'^ = yP^ 



/^/^ 
V 



J 



(A5) 



and a four-momentum integration reads: 

00 1 



11 

" (2^ 2 

p 0-1 



dz \pV 



1 2ir 

dy J dcj). 

1 



We frequently refer to quantities that are 'transverse' 
with respect to a four-momentum Q. The precise mean- 
ing of that is expressed by the transverse projector 



(A6) 



where the normalized momentum is Q'* = I \/ Q^j so 
that 7— matrices that are transverse to Q read: 



7^=T^V-7^-0^'^. 



(A7) 



They satisfy {7^,7^} — 2Tq and 7^7^ = 3. Similarly, 
momenta transverse to Q are defined by 



Pt 



P 



(P-Q)' 



(A8) 



Appendix B: Current diagrams in the three-quark 
framework 



In this appendix we detail the construction of the nu- 
cleon's axialvector and pseudoscalar current in rainbow- 
ladder truncation which is illustrated in Figs. 2 and 7. 
We first collect the properties of the nucleon's covariant 
bound-state amplitude that enters the current diagrams. 
The nucleon amplitude including its full Dirac, flavor and 
color dependence reads 



\n=l 



£abc 



(Bl) 



where \^ n\ai}-i&{Pi P) is the spin- momentum amplitude 
and [Fnjabcd are the flavor tensors. Both of them trans- 
form as doublets under the permutation group (hence 
the indices n — 1, 2), such that the bracket in Eq. (Bl) is 
a permutation-group singlet. The Dirac amplitudes 
carry three spinor indices a, 7 for the quark legs and 
one spinor index S for the nucleon, and they are mixed- 
antisymmetric (5'i) or mixed-symmetric (^"2) under ex- 
change of the indices a, j3 and related quark momenta. 
They can be decomposed in a basis of 64 orthonormal 
Dirac structures X^^^^ which correspond to s— , p— and 
d— waves in the nucleon's rest frame: 



['^n]aPl5{p,q,P) 



64 

Y.myQ^p^s{p,q.P). (B2) 

i=l 



see Ref. [18] for details. The Lorentz-invariant dressing 
functions {t) are the solutions of the covariant Faddeev 
equation, and t abbreviates the five Lorentz-invariant 
momentum variables that can be constructed from the 
momenta p, q and P {P^ — — is fixed). The dressing 
functions for n — 1,2 are not independent but related via 
permutation-group symmetry. 

Similarly, the two isospin-i/2 fiavor tensors F„ carry 
three isospin indices a, b, c for the quarks and one for the 
nucleon. They are mixed-antisymmetric (Fi) or mixed- 
symmetric (F2) with respect to a, b: 



[Fl]abcd = ^ [»a'2]ab (^cd , 
[F2]abcd ~ *f2]ab CTcd , 



(B3) 



where the ak are the Pauli matrices, and they are nor- 
malized to unity: [Fjj,]bad'c[F„]abcd = (^n'K (^d'd- To projcct 

onto proton or neutron flavor states, the nucleon index d 
must be contracted with either of the two isospin vectors 
(1,0) or (0,1), respectively. Finally, the antisymmetric 
tensor sabc in Eq- (Bl) denotes the color part which we 
also normalize to 1. 

The rainbow-ladder truncated current is the sum of the 
impulse-approximation and kernel diagrams of Fig. 7. In 
principle one must also add up all three permutations 
a = 1,2,3, where a is the label of the quark that cou- 
ples to the current. However, the symmetry properties 
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of the nucleon amplitude relate the permuted form-factor 
diagrams among each other via Eqs. (B13) below, such 
that it is sufficient to consider only one of the three per- 
mutations explicitly. In the following we choose a = 3, 
which corresponds to a coupling to the upper quark line 
as shown in Fig. 7. We discuss the construction of the ax- 
ialvector current J^; the equations for the pseudoscalar 
current J5 are completely analogous. 

Starting with the flavor traces, the axialvector and 
pseudoscalar qq vertices come with flavor factors (7^/2 
whereas the rainbow-ladder kernel is flavor- independent. 
Hence, the flavor trace in both a = 3 diagrams is given 
by 



p(3) 



k 

d'd 



it 



bad'( 



Y 



abed ■ 



(B4) 



We keep the doublet indices n' and n for the outgoing 
and incoming nucleon amplitude general so we can treat 
the Dirac and flavor parts in the current separately. The 
above expression vanishes if n ^ n' , and for 11 = n' it 
yields: 



r(3) 
11 



l(3) 
22 



(B5) 



Next, the color trace in the impulse-approximation dia- 
gram equals 1 while the kernel diagram picks up a color 
factor —2/3. The current for the coupling to quark a = 3, 
which is a matrix in Dirac and flavor space, can then be 
written as 



t(3) ■ 



J. 



(3), IMP 



2 t(3),K 

3 ''n'n 





p(3) " 


S'S 


n'n 



d'd 



(B6) 



where the first bracket includes the Dirac parts. We sup- 
pressed the Dirac and isospin indices on the left-hand side 
and the index /i only appears in the axialvector case. The 
total current of Eq. (45) is then the sum over the three 
permutations and the final and initial doublet configura- 
tions: 



a—1 n'n 



(B7) 



Let us now examine the Dirac contributions in Eq. (B6) 
in more detail. The impulse-approximation diagram 



reads explicitly: 



" t(3),IMP' 


Irll 







n' a' S''-y 



P Q 



{pi,qt,Pi) , 

(B8) 

where pi, qi and p/, g/ are the incoming and outgoing rel- 
ative momenta; Pi and Pf are the incoming and outgoing 
nucleon momenta; pi, p2, Pz and p^ = pz±Q/2 are the 
quark momenta; S{pk) are the dressed-quark propaga- 
tors; is the dressed axialvector (or, by analogy, pseu- 
doscalar) vertex; and p and q are the two loop momenta. 
For a symmetric momentum-partitioning parameter 1/3 
the relative momenta are explicitly given by 



Pf =p 



Q 
3 ' 



Pi=P 



3 ' 



and the quark momenta by 



Pi 



P 
2 



P 
"3 



P2 = q 



<if 



p 
2 



P 
3" 



(B9) 



(BIO) 



r(3),K 



and Ps = p + P/3. 

The second diagram involving the rainbow-ladder ker- 
nel is identical to Eq. (B8) except that the incoming am- 

(3) 

plitude \l/„ in the third row is replaced by : 

1 = [Eq.(B8)]^ ^^(3), (Bll) 

do n n 

and the latter is defined by 

[^n^]al3^siPi,qt,Pi) = J Kaa'f3P'{k) X 

(B12) 

X Sa'a"{Pl - k) Sp'p"{p2 +k)x 

{pi.qi + k,Pi) . 

Since Eq. (B12) is one of the three diagrams in the Fad- 
deev equation and has the same structure as Eq. (Bl), it 
is actually not necessary to compute it anew in the form- 
factor calculation. Using the labels / and i to abbrevi- 
ate the final and initial amplitudes, the Dirac bracket in 
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Eq. (B6) can be schematically written as 



2^(3) 
3 



(3) 



[^3r 



2 (f> 

3 



(3) 



Putting the pieces together and adding up Eq. (B13) to 
obtain (B7) yields the following result for the total isovec- 
tor current: 



S'S 



2 



d'd 



(B14) 



(3) 

where Jnn are the Dirac contributions to the a — 3 dia- 
grams, i.e.: 



where $i is the wave function S1S2S3 ^'i, and $ 

,(3) 



(3) 



is the 



wave function obtained from "ifl Both of them can 
be collected beforehand when solving the Faddeev equa- 
tion and implemented in the form factor diagram sim- 
ply through evaluation at the proper incoming momenta. 
Then, only the product Eg, which is independent of 
the loop momentum q, needs to be computed explicitly. 

Finally, in order to obtain the full current in Eq. (B7), 
we need to include the permuted diagrams for a — 1,2. 
They can be inferred from the a — S diagram by applying 
the two-dimensional matrix representations A4' , A4" of 
the permutation group S'^, 



M' 



If -I -V3 
2 ^ ^3 -1 



M" 



-1 

-1 



which act upon the doublet indices n', n via [18] 



''n'n ^ 
t(2) _ 



(B13) 



7(3) _ t(3),IMP _ 2 t(3),K 
'Jnn ~ '^nn 3 "^rin 



(B15) 



In the final step we want to extract the form factors 
Gl and G^ that appear in Eqs. (45-46) from the 
axialvector current in Eq. (B14) and its pseudoscalar 
analogue J5. The flavor factors CTfe/2 on both sides of 
the equations are the same and can be factored out. The 
form factors are then obtained from the Dirac traces 



G^ = ^^(^TrK757^}, 
G. = ^Tr{«7a}, 



(B16) 



G5 = ^Tr{J5 75}, 



where Q is the normalized momentum transfer, 7^ is 
transverse to Q, and t = (5^/(4Af^). 
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